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Abstract 

Oscillations in quantum phase about a mean value of vr, observed across mi- 
cropores connecting two '^He — B baths, are explained in a Ginzburg-Landau 
phenomenology. The dynamics arises from the Josephson phase relation, the 
interbath continuity equation, and helium boundary conditions. The pores 
are shown to act as Josephson tunnel junctions, and the dynamic variables 
are the inter bath phase difference and fractional difference in superfluid den- 
sity at micropores. The system maps onto a non-rigid, momentum-shortened 
pendulum, with inverted-orientation oscillations about a vertical tilt angle 

(j) = TT, and other modes are predicted. 
PACS: 03.75Fi, 67.57.-z,74.50.+r 
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The search for Josephson-hke weak-hnk effects in superfluids has a long history |T]- 
In order to see Josephson-hke phenomena, the characteristic length-scale of the weak-link 
should be comparable to the temperature-dependent coherence length C,{T). Thus a can- 
didate superfluid is ^He — B , with a relatively large zero-temperature coherence length 
^(0) ~ 65nm In fact the oscillatory displacement of flexible membranes walling off 
■^He — B baths connected by micropores, induces a Josephson-like, periodic current-phase 
relation Recently, Davis, Packard, and collaborators have observed a remarkable phe- 

nomenon: metastable oscillations of ^He — B with an average phase- difference of vr across 
the weak link 0. The explanation of these, and other complex tunneling oscillations is of 
considerable general interest, since they provide examples of novel macroscopic quantum 
effects. 

Bose-Einstein condensate (BEC) tunneling of neutral atoms between double-well traps 
has been predicted M to support vr-state oscillations, and other modes. The coupling be- 



tween the Ni^2 condensate atoms in wells 1,2, with phase difference (f) is — v^ViiN^cos^ ~ 
— (1 — z^)^/^cos0 where z = (N1 — N2)/ {N1 + N2) is the fractional population imbalance. The 
coupling energy corresponds to [|IU| a non-rigid pendulum PJTT[] of tilt angle 0, momentum 
Pff) = z and length (1 — p^^)^^^. The non-rigid pendulum is shortest when moving fastest, 
and can thus support inverted-orientation oscillations about (p = tt (concave downwards). 
The BEC mapping onto a non-rigid pendulum dynamics is the generic consequence of wave 
function phase factors and the two-state nature of the tunneling: properties that it shares 
with the ^He — B system. However, the BEC dynamics cannot be naively taken over ||12|| : 
features of the constrained helium wavefunction in a pore geometry, must play a role. 

In this Letter we show that the helium wave function boundary conditions in a pore 
geometry, in conjunction with the Josephson phase relation and the continuity equation, 
determine the ^He — B tunneling dynamics, that is indeed similar to non-rigid pendulum 
dynamics, but with a (momentum-dependent) torsion-bar [jl3| due to hydrodynamic flows 
outside the pore Thus dynamical vr-states of temperature-dependent frequency can 
occur, with a rich variety of other modes. 



The wavefunction is depressed by the boundary conditions {be ) inside, and just outside 
a pore length, that is shown to act as a tunnehng barrier, for sufficiently narrow pores. 
Displacements of the flexible membrane wall are proportional to a fractional shift z of the 
depressed wavefunction, that acts as an oscillating 'Josephson piston': it conveys pairs 
through the tunneling (~ i) determined by the (sine of) the phase difference 

0j across the pore. We now derive the (quasiequilibrium) ^He — B dynamics. 

Free Energy: The ^He — B order parameter [§] can be written as $^_j(r) = 
\E'(r)i?^^j(n(r*)),where \E'(r*) = | \E' (r ) | e*'^'-'^ and |\l'(r)| is proportional to the gap, varying 
on a coherence length ^(T). Here B(ri(r')) is a tensorial factor,with 'ri(r)' representing 
axis rot at ions, varying over a textural healing length Itex,^ 5fim PJ.The Ginzburg-Landau 
(GL) free energy is a scalar, with traces over the tensorial products: -F({$(r), V$(r)}) = 
Fo({$(r)}) + -Fgrad({V$(r)}) p. Unlike the superconductor be of vanishing normal gradi- 
ents of the gap, we must use here 'helium' be , with the order parameter vanishing at the 



walls, . This be , supported by experiment, has been used elsewhere JT^ ,and is 

particularly relevant in the system regime |^-|^ of pore scales ~ ^, ^ Itex- The be is enforced 
by \E' over length scales ~ ^, and the contributions to Fg^ad from |V\E'| contributions 



will dominate those from |V-B| ~ |Vn| ITH]. The n(r*) textural variables can be set equal 



to their fixed,bath value, f^. Thus Fgrad — ^ -^c/ra(i({V\I'(r)}), and the effective free energy 
describing the pore region is 



d^r 



|V$(rOr + ae|$|^ + -|^| 
2m* 1 



(1) 



where a, 6, m* are GL coefficients (constants in space and time) absorbing tensorial-factor 
traces of order unity; and e = (T/Tc — 1), where is the transition temperature. 

The experimental geometry is a superfluid of total mass density p in a pillbox of 
cross-section A, and length L, closed at one end by a flexible membrane (of spring constant 
C) and at the other end by a perforated wall with Nj = 65 x 65 pores of geometric cross- 
section AjQ and length Ljq. The pillbox bath is immersed in a larger superfluid bath. See 
Fig l.We consider a single pore. The x coordinate is along the pore axis, and the system 



cross-section A{x) varies between A and Ajq. 

To obtain an effectively ID model, we average over the transverse cross-section of the 
system, and scale in the equilibrium wave function |\E'o| =('^kl/^)^''^ where the effective 
superfluid mass density is Ps(T) = m*|\E'oP ~ Thus the ID gap wavefunction is ip^x) = 
[^(r)]/|^o|,where [...] = J d\{...)/A{x). With ^{Tf = / {2m* a\e\) , the free energy from 
Eq.(l) is: 

1, ^, 



F ~ eo y dxA{x) 



(2) 



where the energy density is eo = I^E^opc^lel ~ The averaged transverse gradients U{x) = 
[\dy'^{f)\'^+\dz'^{r)\'^]/\ [^(r)]p are ~ ^'^/Ajo in the pore region, and negligible, ~ ^"^/A in the 
bulk. Thus U{x) is an effective energy barrier, arising from the be . We set U{x) = 
in an effective barrier length |x| < Lj/2, and zero outside it. The barrier length Lj is greater 
than the geometric pore length Ljo, as the wave function depression (and large transverse 
gradients) will persist by continuity, for a small distance or 'overhang', outside the pore 
openings, of size 6 = (Lj — Ljo)/2 as shown in Fig 1. We treat 7, > 1, as a fitting parameter, 
and consider e~^^' <^ 1 as small. 

Wavef unctions: The barrier region solutions of the (linearized) GL equation derived from 
Eq. (H) are e^'^^, with a real decay wavevector K = {^"^ /Ajq) — ^"2, for ^/Aro < 7^.( An 
estimate, for square pores, is 7 = v^tt.) With a phase difference 0j = 0(x = Lj/2) — 
(f){x = —Lj/2) across the junction, the (dimensionless) pore wavefunctions are in a form |^ 
generalized to 

^(x) = [a^e-'^'/'^smhK{{Lj/2) - x) + 

a+e^'^^/^sinh K{{Lj/2) + x)]/sinh KLj, (3) 

where coefficients a± here match on to bath wave functions, rather than plane waves [Q. 
Substitution into Eq. (E[) yields a leading-order, |a;| < Lj/2 contribution 



Fj = -EjV(aV-)/«oCos0j, (4) 
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neglecting O {e~'^^^' cos 20 j) 

corrections. The Josephson energy is Ej = 2AjQK^^eoalf(T)/siiah.KLj, and /(T) is a 
possible renormalization factor discussed later. We now estimate a± amplitudes and their 
dependence on interbath pair transfers. 

Equilibrium wave-functions deep in the baths 1,2 are rigid, as bulk ^He — B is incom- 
pressible. Solutions of the nonlinear GL equation from Eq.(l) are hyperbolic tangents, 
fiat in the bulk and falling to zero with slope 1/V2^ at walls fl]]. We take the averaged 
bath wavefunction as 'il>i,2ix) = e^'^'^^/'^g{^x — Ljo/2), where g{x) = tanh(x/v^^). At 
X = ^Lj/2, continuity with Eq. (|^) yields equal equilibrium pore-amplitudes a_ = a+ = oq, 
and = go, where the bath wavefunction at the barrier edge is go = g{S). See Fig. 1. 
Slope-matching determines the (weakly T-dependent) overhang S as sinh(-\/25/0 = V^/K^ 
coth.{{K{Ljo+26)/2). For long narrow pores, go ~ S/ V^^, and the overhang 5 ~ 1/K ^ A^J^ 
vanishes, as pores close, Ajo — > 0. 

Slowly oscillating and small {X{t) ^ ^) displacements of the far-off flexible membrane 
will be transmitted by the rigid bulk, to induce shifts by ri{(x X) of the wavefunction near 
the pores. Since r] will be related to the tunneling through the pores, we term the movable 
region of length Lj in Fig. 1, a 'Josephson piston'. The bulk wavefunctions at x = ^Lj/2+r] 
then match onto the (shifted) piston Eq. (§), yielding a±. For small r], defining z = 27]/ 5 as 
a fractional piston-overhang shift, we can expand as \Ja\a?_/a^ ~ \/l — z'^ to obtain Eq. (^ 
as Fj = Fj{(f), z). (Additional wavefunction stiffness terms, Fj ~ +z'^ make z = the piston 
rest position; these are dominated by membrane stiffness terms, given later). 

The shift of the interfaces between the Josephson region and the bath to x = —^Lj + 
r], X = \Lj + 1] means that there are increases/decreases in the number of atoms in the 
two tunneling overhang regions. Note that z can be interpreted as a "BEC-like" fractional 
population imbalance |p of the overhang volume 2Ajo5, and is not a fraction of the total 
Cooper pairs |p!2| . 

Clearly, the incremental number change (ini,2in each bath will be the number density 
at the interface times the incremental shift drj. The change in tunneling population dn = 

5 



^{drii — dn2) from the overhangs on the two sides is then [0] dn = —{pAjqQqS /2m^)d. 



iz 



where is the hehum atom mass. The piston displacements induce a reactive velocity 
component |JT6|] of the membrane X oc |((7g5)i proportional to the tunneling current. 



Current Equation:The tunneling current (piston velocity) equation for i is provided by 
the continuity equation p(r) = V- J(r) where J is the total mass current density. Integrating 
over baths 1, 2 (outside the piston regions) and using Gauss's theorem, m^h = ^Ajq{[Jx{x = 
—Lj/2)] + [Jx{x = Lj/2)]). Assuming no normal fluid flows through the narrow pores, and 
with superfluid mass-current density [Ja;(a;)] = ps:^Irn{ip*dip/dx), Eq.(|) yields: 

— V 1 — -2 sm0j = — — — . (oj 



2Kjdt d(t)j 
where H and Kj are defined below. The physical picture is of phases at Nj pore openings 
in each bath locked together, and with a common difference between baths,so that the 4000 
Josephson pistons vibrate in unison. 

Phase Equation: The Josephson relation for = 02 ~ 0i) the total interbath phase 
difference, is = A/i//i and the chemical potential difference from the tunneling transfer of 
n atoms will contribute as A/i ^ —dFj/dn.{T]ie bulk contributions from the incompressible, 
constant-density fluid in the baths, clearly cancel out of the difference A/i.) Work done by a 
fixed [0 external pressure Pext, and by the Josephson piston, pushing against the membrane, 
of spring constant C, can also be included |16[. Thus 



h dcf) Pext z dH 

Az H — - — I — -j^=^= COS0J = — — ; (baj 



2Kjdt Pj ^/Y^ ^ dz 

A 2 I Pext^ 



H = -Az^ + - cos 0j, (6b) 

2 Pj 

where Pj = pKj/m^ is a pressure scale. The frequency scale is 2K,j/h = 
ipJp)i2hK/m*5)fiT)/smHKLj) ~ |e|/(T); and A = (C g^d) / i2APj) ~ C//(T) is a (di- 
mensionless) membrane stiffness parameter. 

Note that and z are not canonically conjugate, as 7^ 0j due to hydrodynamic effects 
The hydrodynamic bath-pore superfiow / = {Ah/ AjQ)AjQPs{h/m*){(j)h/Lh) is driven 
by a phase gradient (ph/Lh over a hydrodynamic length Lh in each bath. (Here Ah/Ajo is 
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a bath-pore 'transmission probability'). The same current / = /jv^l — sin tunnels 
through the pore, where the Josephson critical (mass) current is Ij = Ej2m'i/h ~ |ep/(T). 
With total phase change = 2(j)h + 0j, we have a nonlinear relation |^ generalized to 



+ aVl — sin0j, (7) 



where a = {2'^){KLhgQf{T))/smh.{KLj) ~ |e|/(T). This corresponds to a torsion-bar as 



in the inductive SJJ case |110| , p!3| , but now, momentum-dependent. 

The temperature-dependence of A(T) is essentially that of /(T). Experimentally,the 
critical current goes as ~ (1 — T/Tc^)^, appearing to vanish before transition, at a Tea = 
0.91Tc §. Since Ij ~ e'^f{T), this 'experimental' form /(T) = (1 - T/T^af/e^, could be 



used. Or, /(T) could be attributed to thermal phase fluctuations fl^.We estimate it, in this 
picture [jl8[. 

Eq. (^,^J3) constitute the model equations for the ^He — B system, in terms of the mu- 
tually linked dynamics of the fractional Josephson piston shift t]/6 = z/2 and the phase 
difference 0j across it. They can be solved in terms of elliptic functions [|r^] just as in the 
a = 'BEC-like' case [^.They correspond to a non-rigid pendulum of tilt angle 0j,and 
angular momentum 2;, with a momentum- dependent torsion bar. The torsion-bar [l^ in- 



duces hysteretic effects for a > l.We focus on a < 1 where the modes are essentially as for 
a = , with the equations predicting five distinct undriven modes (instead of the two of the 



rigid pendulum). The modes can be visualized by plotting [|rT| the locus of the pendulum 



coordinates {^/l — z'^ sin 0j, — — z"^ cos 0j). With a dc drive Pext there is an ac Joseph- 
son oscillation at a frequency Uac = 2m3Pext/ P^', with added resonant ac drives, there are 
Shapiro-like dc current spikes P,p!6|. 'Zero-state' oscillations of the pendulum can occur, 



with time averaged (0) = = (z).The pendulum, once excited, can rotate freely, with a 
running phase, and a 'self-trapped' piston-shift /membrane-displacement ~ (z) 7^ 0. 

Since the non-rigid pendulum is shortest, when moving fastest, it can execute small 'in- 
verted' oscillations about an average value of (0) = tt. These A < 1, 'vr-state oscilla- 



tions', with {z) = 0, and are dynamically metastable, with the momentum-shortened pen- 
dulum 'digging a well for itself by its motion. Finally, there are two types of A > 1, 
'vr-state rotations' with self-trapped (z) ^ and (0) = 7r,(on either side of a fixed point 



(j) = 71, z = Vl — A^2). These correspond to the non-rigid pendulum executing closed-loop 
trajectories, floating above the point of support [§,|11|]. Damping of vr-state oscillations 
through the non-rigid pendulum velocity , or the momentum z, behave quite differently. 
While the former drives the system to rest at 0j = 0, the latter damps it to 0j = vr. 



We use values [|] Ljo = 0.5 10'^ cm, Ajo = 10"^° cm^, p = 0.08 g/cm^, m* ~ 2m3 = 
10-23 = 4225, ^(0) = 65 nm, A = 0.7 cm^, C = 10^ dynes/cm, = 0.91 mK, and a 

purely illustrative parameter 7 ~ 23 . Then 6 ^ 90A and 2Kj{0)/h ~ 400 Hz, {lj{0)/ms) ~ 
3.2 X 10^ atoms/sec, corresponding to a Josephson energy per pore Ej{0) ~ 0.02 ksTc. We 
find A(0) ~ 0.9; Pj - 0.35 mPa; and a(0) ~ 0.5. 

The results of our model are shown in Fig.2. The dimensionless membrane displacement 
versus dimensionless time {t2Kj/h t) is shown for zero- and vr-states. The insets show 
the temperature-dependent frequencies and critical currents. These include: i) zero-state 
oscillations , with harmonic frequencies uJoiT) = [(1 + A)/(l + a)]^ {2Kj/h) that vary from 
kHz, to zero at T^. ii) vr-state oscillations, with smaller amplitudes and frequencies uj.,^{T) = 
[(1 — A)/(l — a)Y^'^{2Kj/h) < uJo{T), that vanish on warming, when A(T) crosses unity, 
well before Tc. iii) A critical current [|l^ 20 picogm/s) goes as , appearing to vanish 
before T^. iv) With a phase angle defined as the time-integrated membrane displacement 
0, the current-phase relations (not shown) have same-sign, positive slopes in the zero- and 
vr-states. 

Our model clearly reproduces experimental features of the magnitudes/temperature de- 
pendences,of the oscillations, their frequencies, and currents [§],§[. It also predicts self- 
trapped vr-states with toj^ = [(A'^ — l)/{l — a/A)]^{2Kj/h), that have not yet been reported. 

Further work could include experimental investigations of predicted intermode transitions 
with frequencies dipping to zero, as A(T) or z{0), (f)(0) are varied through critical values 
PJT6[|; and a search for self-trapped vr-states. 
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TT states have been related to a static metastable minimum induced by spin-textures 
in the baths This could dress our dynamic vr states from oscillations of the density 

depression at pores, producing fine structure in the oscillation modes. In fact, bistability 
between two types of tt oscillations has been reported pl[] . 

In conclusion, we have modelled the dynamics of ^He — B baths closed by flexible walls, 
and connected by micropores. The dynamic equations arise from helium boundary condi- 
tions,the continuity equation ,and the Josephson phase relation. The dynamical variables are 
the phase difference, and the shift of the pore wave-function overhang, beyond the geometric 
pore length. The system maps on to a non-rigid momentum-shortened pendulum dynamics. 
The TT— phase oscillations observed in ^He — B have features similar to those predicted by the 
model. Other tunneling modes are predicted. 

It is a pleasure to thank J.C.Davies, R.E.Packard, E.Granato,M. Mehta, R.Simmonds 
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FIGURES 

FIG. 1. Schematic diagrams (not to scale) of : a)Effective pore/bath geometry, showing geo- 
metric pore length Ljo and pore (bath) area Ajo (A). Dashed line shows the rest position of the 
'Josephson piston' of length Lj = Ljq + 25.The piston shift is ij and the membrane displacement 
is X{(x 77). ; b)Equilibrium wavefunction ip{x) versus x. 

FIG. 2. Dimensionless z{t) (proportional to membrane displacement) versus dimensionless time 

for |€| = 1, and parameters as in text, vr-state oscillations (early times); are kicked out to zero 
phase oscillations (later times). Left inset: Zero-state (7r-state) angular frequencies u!o{T) {uJt^{T)) 
defined in text, versus temperature. Right inset: Characteristic current Ic oc Ij{T) the Josephson 
critical current, as defined in text, versus temperature. 
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